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Abstract. The dynamics of the spin-boson Hamiltonian is considered in the stochastic 
approximation. The Hamiltonian describes a two-level system coupled to an environment 
and is widely used in physics, chemistry and the theory of quantum measurement. 

We demonstrate that the method of the stochastic approximation which is a general 
method of consideration of dynamics of an arbitrary system interacting with environment 
is powerful enough to reproduce qualitatively striking results by Leggett at al. found earlier 
for this model. The result include an exact expression of the dynamics in terms of the 
spectral density and show an appearance of two most interesting regimes for the system, 
i.e. pure oscillating and pure damping ones. Correlators describing environment are also 
computed. 
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1. Introduction 



The so-called spin-boson Hamiltonian is widely used in physics and chemistry. In 
its simplest version it describes a dynamical model of a two-level system coupled to an 
environment. One of the bàsic ideas is that the environment induces dissipative effects, 
but as we shall see, the picture is much richer. Examples include the motion of defects 
in some crystalline sòlids, the motion of the magnètic flux trapped in an rf SQUID ring, 
some chemical reactions, some approaches to the theory of quantum measurement and 
many other quoted in the survey paper [1] to which the present work is inspired. The 
"spin-boson" Hamiltonian considered in the present paper is the same considered in [1], 
i.e. 

H x = - l - Aa x + l - ea z + J dku(k)a+(k)a(k) + \a z (A(g*) + A+(g)) (1.1) 

where a x and a z are Pauli matrixes, e and A are real parameters interpreted respectively 
as the energy difference of the states localized in the two wells in absence of tunneling and 
as the matrix element for tunneling between the wells. We set A > and deno te 

A+(g) = J a+(k)g(k)dk , A(g*) = J a(k)g*(k)dk 

where a(k), and a + (k) are bosonic annihilation and creation operators 

[a{k),a+{k')] = ó(k-k') 

which describe the environment. 

We denote u(k) the one-particle energy of the environment and assume uj{k) > 0. The 
function g(k) is a form factor describing the interaction of the system with the environment, 
À is the coupling constant. It is well known that, in times of order í/À 2 , the interaction 
produces effects of order t. Thus À provides a natural time scale for the observable effects 
of the interaction system-environment. 

In the paper [1] it was found a very rich behavior of the dynamics of the Hamiltonian 
(1) ranging from undamped oscillations, to exponential relaxation, to power-law types of 
behavior and to total localization. Leggett et al. [1] found the remarkable result that main 
qualitative features of the system dynamics can be described in terms of the temperature 
(i.e. the initial state of the environment) and of the behavior, for low freqüències u, of the 
spectral function 

J(u) := J dk\g(k)\ 2 5(u(k) - u) (1.2) 

The goal of this paper is to investigate the dynamics of the Hamiltonian (1) in the so called 
stochastic approximation. The overall qualitative picture emerging from this approach is 
similar to the one described in [1] and in some cases also the quantitative agreement is 
good (cf. Section [4]). 
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2. The stochastic approximation 



The bàsic idea of the stochastic approximation is the following. If one has a Hamilto- 
nian of the form 

H X = H + XV (2.1) 
then, by definition the stochastic limit of the evolution operator 

U w (t) = e ltH °e- ltHx (2.2) 

is the following limit (when it exists in the sense specified by (2.11) and (2.12) below): 

Notice that, on the right hand side of (2.3), there is not but its rescaled version 

U^ x \t/\ 2 ). Thus the limiting evolution operator U(t) (2.3) describes the behavior of 
the model in the time scale described in the introduction. The stochastic approximation 
is a natural generalization of the Friedrichs-van Hove limit which uses the same time 
rescaling but allows only to compute vacuum expectation vàlues of the form (U^ 2 Auj;^t) 
for particular classes of observables A. This leads to irreversible evolution and to 
the corresponding master equation. On the contrary, the stochastic approximation 
leads to reversible, unitary evolution and to the corresponding quantum stochastic 
differential equation from which the master equation are deduced by a now Standard 
procedure which consists in integrating away the environment degrees of freedom. 

The stochastic approximation to the original dynamics (2.2) consists in the computa- 
tion of the limit (2.3) in the sense of matrix elements over some states ip\ (called "collective 
states") which themself depend on the parameter À in a singular way. The fact that one 
cannot expect the limit (2.3) to exist for arbitrary states, but only for a carefully chosen 
class of states was already pointed out in the classical paper of van Hove [2] . The effective 
determination of this class of states was obtained in the paper [3]. 

The stochastic approximation could also be considered as a new kind of semiclassical 
approximation in the sense that it studies the fluctuations around the classical solution 
and not the approximation to it. This interpretation however shall not be discussed here 
(cf- [4]). 

One of the important features of the stochastic method is its universality . The re- 
striction to Pauli matrixes in (1.1) is unnecessary: the theory is applicable whenever the 
evolution operator U^(t) (2.2) satisfies the equation 

^^& = -i\V(t)U W (t) (2-4) 
where V(t) = e ltH "Ve~ ltHa has the form 

V(t) = £(£>+ ® A a (t) +D a ® A+(t)) (2.5) 
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and the D a are operators describing the system. The rescaled evolution operator (p-) , 
associated to (2.5), satisfies the equation 

Í^Mi ^(Kv^i*)^*^^))^*) (2.7) 

In the spin-boson Hamiltonian (1.1) the D a are Pauli matrixes (cf. formulae (3.10b), 
(3.12b) (3.8)) and 

A a (t) = í a{k)e- UüJ ^ k) g*{k)dk (2.6a) 



where the functions u> a (k) have the form 

w a (fe) = Lü(k) — ui a ; cü = 1,2,3 (2.66) 

here ui{k) is as in (1.1) and the ui a are characteristic freqüències given by formula (3.12a). 
^From (2.7) it is clear that to have a nontrivial limit for U^ x \t/\ 2 ), the limit 

£3 = (2-8) 

should exist. It can be proved (cf. [4]) that the limit (2.8) exists for "good" functions 
Lü a (k) and g(k) in the sense that 

&o ( X^ (I) ••·1 A - = <£i(*i) ■■■£»(*»)> (2-9) 

where the indices £; label the creators (e = 0) and the annihilators (e = 1); the brackets 
in (2.9) denote mean vàlues over the Fock vacuum or a temperature state; and for each a, 
b a (t) is the Fock Boson quantum field described by (2.14), (2.15) below. In the literature 
5-correlated (in time) quantum fields are often called quantum noises. In the present 
paper only quantum white noises shall appear, but it is important to keep in mind that 
many other possibilities can arise from different physical models. 
^From (2.7) one has in the limit À — > 0: 

=-iY.i D t®^{t) + D a ®bt{t))U{t) (2.10) 

a 

The limit (2.3) means that: 

lim(* A , U™ (j^j ^ = (</,, U(t W ) (2.11) 
where the collective vectors ~^>\ are defined by 



and converge to the corresponding n-particle vectors in the noise space, given by: 



V = 6í 1 (*i)···6í B (*»)V' (0) (2-13) 

and ^ (0) are the vacuum vectors in the corresponding Fock spaces. If u a ^ up 
for a 7^ P (as it is the case for the Hamiltonian (1.1)) then b a , 6+ satisfy the following 
commutation relations 

[b a (t),b+(1?)]=6 aa ,J a 6(t-1?) (2.14) 
where J a is the spectral function (1.2) 

J a = 2n J dk\g(k)\ 2 5(u a (k)) (2.15) 

Thus, as announced in the introduction, in the stochastic limit the spectral function 
emerges naturally as the covariance of the quantum noise. Some care is needed in the 
interpretation of equation (2.10) because, as it is clear from (2.14) the b a (t) are not bona 
fide operators but only operator valued distribution. In order to give a meaning to equa- 
tion (2.10) (more precisely to its matrix elements in the n-particle or coherent vectors), 
we rewrite (2.10) in normal form by bringing b a (t) to the right of U(t). This gives rise to 
a commutator which can be explicitly computed. The result is: 

^ = -i J2(D+U(t)b a (t) + D a b+(t)U(t) - i la D+D a U{t) (2.16) 

a 

where 7 Q are complex numbers given explicitly by: 

7a = f dr í dke lTÜJ ^\g(k)\ 2 (2.17) 

J — oo J 

The connection between the constants 7 Q in the last term in (2.16) ( the Ito correction 
term), and the spectral function (2.15) is obtained by exchanging the <ir-and the dk- 
integral in (2.17) and using the known formula 



r° i 

/ e ltüJ dt = n5{u)-iP.P.- 

J-oo U 



where P.P. denotes the principal part integral. This shows that the spectral functions are 
the real parts of the constants 7 a; emerging in the Ito correction term. This connection is 
the prototype of the dispersion relations widely used in quantum physics since its origins. 
Since the 7 a are complex, equation (2.16) looks like an equation driven by a non self- 
adjoint Hamiltonian. However this is only an apparent phenomenon due to the normal 
order. The true Hamiltonian (2.10), although singular, is formally self-adjoint and this 
gives an intuitive explanation of the unitarity of the solution of (2.16) or, equivalently, of 
(2.10). 

The relations (2.14)-(2.17) define the stochastic approximation to the system (2.4), 
(2.5). The term stochastic is justified by the fact that the distribution equation (2.16) 
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which has a weak meaning in the n-particle vectors, can be interpreted as a quantum 
stochastic differential equation (and, in fact, it is in this form that this equation was first 
derived [3]). The operators b a (t), b~^(t) are called a quantum white noise and the 
additional term in (2.16), arising in (2.16) from normal order is called the drift or the Ito 
correction term. More precisely, in quantum probability one usually writes (2.16) in the 
form 

dU(t) = -i J2(D+dB a (t) + D a dB+(t) - i la D+D a dt)U{t) (2.18) 

a 

where 

/t+dt 
b a (r)dT 

are called stochastic differentials and satisfy the Ito table: 

dB t dB+ = 2-fdt ; dtdB+ = dB t dB t = dB+dB+ = dB+dB t = (2.19) 

The proof of the Ito table (2.19), as well as its rigorous meaning, was first established in 
[5]. This has been subsequently applied to several models in quantum òptics by [3] and 
[4]. Using it, the unitarity of the solution of (2.18) is easily established. 

The advantage of equation (2.16) over the original one (2.4) is that it is in some sense 
completely integrable and one can easily read the physics from it. For example for the 
vacuum expectation value one has the equation: 



d{U(t)) 



-Y,-yaD+D a {u(t)) 



dt 

a 

which gives the damped oscillatory regime (the 7 a are complex number): 

(U(t))=e- rt , T = J2^D+D a 



In the following we shall apply this method to the Hamiltonian (1.1) and, in Section (4), 
we shall compare our result with those of [1] , [6] . 
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3. The stochastic approximation for the "spin— boson" system 

In order to apply the stochastic approximation to the Hamiltonian (1.1), we write 
(1.1) in the form (2.1) where 

H = Hs + H R (3.1) 



The system Hamiltonian Hs is 



H s = ~Aa x + ^ea z (3.2) 



and the reservoir Hamiltonian Hr is 

H R = J dku(k)a+(k)a(k) (3.3) 

The evolution operator U^ x \t) satisfies equation (2.4) where 

V(t) = a z (t)(A(e-^g*) + A+(e^g)) (3.4) 

and 

a z (t) = e itHs a z e- itHs (3.5) 

To bring (3.4) to the form (2.5) let us compute (3.5). The eigenvalues of the Hamiltonian 
(3.2) are 

H s \e± >= À±|e± > (3.6) 

where 

A± = ±^Az/ (3.7) 

\e±>= A (} ) (3-8) 



and 



e . i íe^ 2 



^± = K ±^, ^ = v 1+ UJ (3 - 9) 



Notice, for future use, that: 



1 - y 2 

{e±\a z \e±) = — — % ; (e + |a 2 |e_) = (e_|a 2 |e + ) = \jv 



l + fji T 
Therefore 



tL. DD+ + ?± 

1 + a»Í 1 + M+ 



a z(t) = , g DD+ + ; P + D+D + 1/"^**^ -D + ^- 1 e- í ^ A J D + (3.10o) 
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where 

D = |e+Xe_| (3.106) 
The interaction Hamiltonian (3.4) can now be written in the form (2.5): 

3 

V{t) = "£(D+ ® A(e- <tw V) + h.c.) (3.11) 

where the three spectral freqüències correspond respectively to the down, zero, and up 
transitions of the 2-level system, i.e. 

uj l {k) = üj{k)-vA ; u 2 (k)=u(k) ; u 3 (k) = u(k) + vA (3.12a) 

1 - u 2 1 - u 2 

D 1 = u~ 1 D+ ; D 2 = DD+ + ^D+D ; D 3 = u~ 1 D+ (3.126) 

1 + n_ 1 + 

The corresponding limiting evolution equation therefore has the form (2.16). It is impor- 
tant to note however that the constants (2.15) for a = 2, 3 vanish, i.e. 

J 2 = h = (3.13) 

We shall see that the purely oscillatory regime, first discovered by Leggett et al. [1] 
corresponds to the case when also J\ vanishes. In this sense it can be interpreted as 
an off-resonance regime. In this regime a strange (from the point of view of stochastic 
theory) new phenomenon take place: in í/À 2 -limit the environment disappears (i.e. the 
limit on the right hand side of (2.8) is zero, corresponding to a quantum white noise of 
zero variance). However a remnant of the interaction remains because, after the limit, the 
system evolves with a new hamiltonian, equal to the old one plus a shift term depending 
on the interaction and on the initial state of the field. This is a kind of Cheshire Cot effect. 

(3.13) implies that the operator s b 2 and 63 should be absent in (2.17). However the 
constants 72 and 73 as well as 71 do contribute to (2.16). We denote bi(t) by b(t). Thus 
the operators b{t), b + (t) satisfy 

lb(t),b + (t')]= 7 5(t-t') (3.17) 

with 7 given by (3.15) below and v (in 7) given by (3.9). The limiting evolution equation 
can then be written: 

= Db+(t)U(t) - D+U(t)b(t) - (7 + ia)D+DU(t) - i<pU(t) (3.14) 

where 

7 = z/" 2 7rJ(z/A) , (3.15) 



a = u- 2 (I(-uA) -I(uA)) + 



1 — fi 



2 x 2 
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„ = „-»I(-„A) + g^ 

and we denote 

J(üj) = í dk\g(k)\ 2 S(u(k) - u) ; I(u) = P.P. H du; ' J ^ (3.I6) 
J Jo w — w 

where P.P. means the principal part of the integral. 

In the notations of quantum stochastic equations (3.14) reads 

dU{t) = (DdB+ - D+dB t - (7 + ia)D + D - i(p)U(t) (3.18) 

Notice that all parameters 7, a and ip in the evolution equation (3.14) are expressed in terms 
of the spectral density J(w) (3.16) and parameters A and e of the original Hamiltonian 
(1.1). 




4. Analysis of the stochastic approximation. Zero temperature 

Let us discuss now in more detail the implications of the results of the previous sections 
for the "spin-boson" Hamiltonian. All the information about the model is encoded into the 
constants 7, a and ip and these constants are expressed in terms of the spectral density J(u) 
(3.16) depending on the parameters of the Hamiltonian (e and A) and the temperature 
(not yet introduced up to now). Thus the method of stochastic approximation confirms 
the conclusion of Leggett at al. [1] that the long-time behavior of the model is expressed 
i n terms of the spectral density J(ui). 

Now let us discuss the dynamics of the system in the stochastic approximation. We 
are interested in a pure damping or pure oscillating behavior. 

For the vacuum expectation value we have 

(U(t)) = e~ l(pt + e -^( e -(7+-)í _ i)D+D (4.1) 

and taking trace over the spin variables one gets (since TrD + D = 1 ): 

(trU(t)) = e-fr+^+v)]* (4.2) 

Since 7, cr, and <p are real (cf. (3.14)-(3.16)), one has a purely oscillating behavior (4.2) if 
and only if there is no damping, i.e. 

7 = (4.4) 
However one cannot have a vanishing of oscillations, because the quantity 

a + p = v- 2 I(-vA)+^ r ^^J 7(0) >0 (4.3) 
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is strictly positive for positive J(uj) (because v, A > 0, cf. (3.9)) and I(ui) is given by 
(3.16). 

The stochastic approximation to the vacuum expectation value of the Heisenberg 
evolution of a z is given by 

P(t) = (U*(t)or z (t)U(t)) (4.4) 
^From equation (3.14) one gets the Langevin equation for P{t) which solution is 

P(t) = v^e-^iD+e^-^* + De- l(a - uA)t ) + 

+D+D - + í-^ (4 5) 

Let us discuss separately the simplest case e = 0. 

The case e = 0, zero temperature 
In this case one has 

P(t) = e-^(D+e l ^- A ^ + De- ï{a - A)t ) (4.6) 
where 7, a, and I(uj) are now 

7 = ttJ(A) ; a = I(-A) - I(A) ; = P.P. I ^1^1 

J U) — U) 



Two interesting regimes can now appear: 
(i) No oscillations . In this case 

a - A = (4.7) 
Equation (4.7) is equivalent to the integral equation 



/ 



- P.P. zní^l = A (4.8) 
x + A x-A 

If equation (4.8) is satisfied then we have pure damping: 

P(t) =e~' ít (D + + D) (4.9) 

We will discuss solutions of eq. (4.8) later. 

Another regime is 
(ii) Pure oscillations . This regime is defined by the condition 

7 = 7iJ(A) = (4.10) 

Notice that, because of (3.16) this condition defines an off— resonance condition. 
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If equation (4.10) satisfied, then 

P(t) = D + e i{a - A)t + De- l{a ~ A)t (4.11) 

where 

a _ A=r í dxJ(x) _pp f dxJ{x) _ A 
a J x + A ' J x- A 

This case of pure oscillations is very interesting. If there is a damping then after a rather 
short time P(t) becomes a small quantity which is difficult to observe. The case of perma- 
nent oscillations looks more promising for observations. This regime is of primary interest 
in the context of the so-called macroscopic quantum coherence phenomenon [7]. 

The purely oscillatory regime was discovered in [1] but the region of parameters there 
is different from ours. To get pure oscillations we need the only off-resonance condition 
(4.10), i.e. in terms of the spectral density what we need is: 

J(A) = J dk\g(k)\ 2 5(u(k) - A) = 

The difference with [1] can be attributed to the different boundary conditions on correla- 
tors. 

Let us present our results on the computation of the correlator 
C{t) = \({Uta z U u a z })= 1 -{P{t),P{Q)} 

We have 

C(t) = I e -(7+-+-A)t ( u -2 + u -i D ( lz]£_ + 

+h . c . + Izg {V - HD + D+) + izg. DD+ + + 

The trace of C(t) is 

(2 2 \ / 2 \ 

— ^± ^ ) ( — ^ ) 

The qualitative behavior of C(t) is such as for P(t). 

Non-zero temperature 

For a non-zero temperature we get a stochastic evolution equation of the same form 
as before (3.14) only with new constants 7, a and ep. More precisely: 

7 = z/- 2 tt(J + (z/A) + J_(i/A)) 
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(4.12) 



i-/4 
i + /4 



1-/J.Í 



(I + (0)+/_(0))+ 



+z/- 2 (/+(-z/A) - J+(z/A) + 7_(-z/A) - I-(uA)) 
where spectral densities are 



J(ü/) 



1 - e -^ w 



j_M = 



J(u)e 



-fiu 



1 - e-#" 



Here J(w) is the spectral density (3.16) and /? is the inverse temperature. 
The functions I±(üj) are defined by 



Í±M = p.p. j _± 



du/J±(u/) 
cu 



One has the same as for the zero-temperature expressions (4.15) and (4.12) for P(í) and 
C(í) but now with new constants 7 and a depending on temperature: 



7 = v 2 7í J(uA) coth 



12 



Conclusion 

To conclude the following main result are obtained: 

The theoretical role of the spectral function is explained through its emergence from 
a canonical limit procedure. Moreover this function is shown to be real part of a complex 
function whose imaginary part defines an energy shift in the system Hamiltonian. When 
the environment free energy depends only on the modulus of momentum (oj(k) = uj(\k\) 
in (1.1)) the real and imaginary part of this function are related by a Hilbert transform, 
thus making a bridge with the Standard dispersion relations (cf. Section (2.)). 

In the stochastic approach not only the Heisenberg equation of the system observ- 
ables is controlled, but also the environment evolution. It is shown that the environment 
converges to a quantum noise (a master field, in the particle physicists terminology) . This 
gives a theoretical (i.e. based on a microscopic Hamiltonian description) foundation to the 
use of classical of quantum noises widely used in several contemporary approaches to quan- 
tum measurement theory [1] , [7] . We can compute the limit matrix elements of Heisenberg 
evolution for arbitrary n-particle or coherent vectors. The vacuum matrix elements give 
rise to the master equation. The control of the other matrix elements is a new feature of 
the stochastic approach. 

The purely oscillatory regime, discovered by Leggett et al. [1] is related here to a 
Cheshire Cal effect in which the environment variables vanish in the limit but the inter- 
action leaves a track in the system behavior in the form of an operator shift in the system 
Hamiltonian (cf. Section (4.)). 
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